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ABSTRACT 

We study the clustering of the mass distribution in cold dark matter models using 
large cosmological N-body simulations. We investigate spatially-flat models with a 
cosmological constant and scale- invariant {n = 1) primordial power spectra, as well 
as open-bubble inflation models. All the models we consider are normalized according 
to the fluctuation amplitude measured in the COB£^-DMR microwave background 
anisotropy data. With an age of the universe to ~ 14 Gyr (12 Gyr) for the flat 
(open) models, a baryon mass density parameter fie = 0.0125/i^^, and a reasonable 
assessment of the systematic uncertainties in the cumulative cluster mass function, 
the observed abundance of rich galaxy clusters leads to tight constraints on the mass 
density parameter JIq- The allowable ranges are 0.4 ^ f2o ^ 0.5 for open models and 
0.25 ^ r2o ^ 0.4 for flat models. The upper limits on fio can be relaxed if one lowers 
the Hubble pararameter and increases the age of the universe, but h < 0.25 is required 
for fio = 1 to be allowed. The constraints also change if one allows tilted primordial 
power spectra. An f2o = 1 cold dark matter model with h = 0.5 can be constructed 
to satisfy both the cluster and DMR constraints, but it requires a tilted primordial 
power spectrum, with n w 0.8 and a corresponding contribution to the DMR signal 
from gravitational waves that reduces the implied cts by a further 27%. We compare 
the evolved mass correlation functions and power spectra of the most promising of 
our N-body models with those of galaxies in the APM survey. The flat models have 
steep correlation functions at small scales and require the galaxy distribution to be 
antibiased on scales r <^ 8h~^ Mpc. The open models require little or no antibias on 
small scales and a positive bias on large scales; these biases are small for il 2± 0.4, 
implying that, in this case, galaxies approximately trace the mass over a wide range of 
scales. The lack of a positive bias on small scales in almost all of these N-body models 
is difficult to reconcile with the mean mass-to-light ratio of cluster galaxies which, 
if r^o ^ 0.2, implies that galaxies are overabundant in clusters relative to the field. 
The tilted fig = 1 model, on the other hand, does require that galaxies be positively 
biased on small scales, and that the bias to become stronger on larger scales. We also 
compute the topology of isodensity contours in these models, obtaining theoretical 
predictions that are less sensitive to the details of galaxy formation. 
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1 INTRODUCTION 

Earlier this century, Einstein and de Sitter (1932) sug- 
gested that until there was sufficient observational data, 
progress in constraining the cosmological model could per- 
haps be best made by focussing discussion on the restricted 



model that now bears their names. In this Einstein-de Sit- 
ter model, the spatial hypersurfaces are flat, the cosmolog- 
ical constant A vanishes, and the universe is dominated by 
pressureless matter, so the mass density parameter flo — 
87rGp6(to)/(3-ffo^) ~ 1, where pb{to) is the mean mass den- 
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sity today, Hq is the Hubble parameter, and G is the New- 
tonian gravitational constant. 

In the 1980s, the ideas of inflation and pressureless, 
cold dark matter led to a more encompassing and far more 
predictive version of the Einstein-de Sitter model: the cold 
dark matter (CDM) scenario (Peebles 1982; Blumenthal et 
al. 1984; Davis et al. 1985). In this picture, inflation ex- 
plains why the universe is homogeneous on large scales, and 
quantum fluctuations generated during inflation provide a 
source of primordial perturbations that can grow through 
gravitational instability of the cold matter into the galax- 
ies, clusters, and superclusters that are observed today. The 
discovery of large-scale anisotropies in the cosmic microwave 
background (CMB) by the COBfi-DMR experiment (Smoot 
et al. 1992) gave a big boost to the broad picture of struc- 
ture formation by gravitational instability since it revealed 
inhomogeneities in the early universe at roughly the level 
predicted by theoretical models. However, the DMR results 
also led to a powerful argument against the simplest version 
of the cold dark matter model: when normalized to produce 
the observed level of microwave background anisotropies, the 
model predicted excessively massive galaxy clusters (e.g., 
Wright et al. 1992; Efstathiou, Bond, & White 1992a; Bah- 
call & Cen 1993; Bartlett & Silk 1993; White, Efstathiou & 
Frenk 1993a). The level of discrepancy is only a factor of two 
or three in the fluctuation amplitude, but it reinforces years 
of accumulating evidence for a difference between the shape 
of the matter power spectrum predicted by standard CDM 
and the shapes of the observed power spectra of galaxies 
and galaxy clusters (e.g., Bahcall & Soneira 1983; Klypin & 
Kopylov 1983; Maddox et al. 1990; Efstathiou, Sutherland 
& Maddox 1990; Eke et al. 1996a). 

In the wake of the COBE results, many variations on 
the cold dark matter theme have been explored. Arguably, 
the most attractive of these models are those that retain 
the assumptions of cold dark matter and an inflationary ori- 
gin of perturbations but drop the assumption that ilo = 1. 
These low density CDM models come in two varieties, spa- 
tially flat models with a cosmological constant (e.g. Efs- 
tathiou, Sutherland & Maddox 1990, Ostriker & Steinhardt 
1995; Ganga, Ratra, & Sugiyama 1996c; Liddle et al. 1996; 
Lineweaver et al. 1997; Peacock 1996; Ratra et al. 1997), 
and open models with A = (e.g., Kamionkowski et al. 
1994; Gorski et al. 1995,1996b; Ganga et al. 1996c; Peacock 
1996; Ratra et al. 1997). One can consider models with both 
space curvature and non-zero A, but the observations have 
not yet driven many theorists to such extremes. Within some 
rather broad ranges of parameters, both flat and open mod- 
els can accommodate the DMR data on their own. A key 
element in tightening the range of acceptable parameters 
has again been the comparison between the DMR normal- 
ization and the predicted masses (or abundances) of galaxy 
clusters. Most implementations of this approach have used 
semi-analytic computations based on the Press-Schechter 
(1974) formalism (e.g.. White et al. 1993a; Borgani et al. 
1996; Eke, Cole & Frenk 1996b; Gorski et al. 1996b; Ki- 
tayama & Suto 1996; Viana & Liddle 1996). In this paper we 
use large cosmological N-body simulations to study the clus- 
ter mass function in DMR-normalized, open and flat CDM 
models with Q,o = 0.1, 0.2, 0.3, 0.4, 0.5, and 1. We also 
compute clustering properties of the non-linear mass distri- 



bution in the most promising of these models and in a tilted 
(n < 1) Jlo = 1 model. 

Variants of the CDM scenario that retain the assump- 
tion of f^o = 1 are still under discussion. These variants 
include the tilted model which we consider, but also may 
involve a broken primeval fluctuation spectrum, a gravita- 
tional wave contribution to the CMB anisotropies observed 
by COBE-DMK, an admixture of hot dark matter, a high 
baryon fraction, a decaying particle that boosts the neutrino 
background, or a low Hubble parameter — in short, one can 
fiddle with inflation, flddle with the matter content, or flddle 
with Ho (see, e.g.. Bond & Efstathiou 1991; Schaefer & Shafi 
1992; White et al. 1996). These models can, to varying de- 
grees, account for the COBE-DMH observations and large 
scale structure, provided that optically bright galaxies are 
biased with respect to mass by a factor 6 ~ 1.5 — 2. However, 
there are two other arguments that generically favour low 
density models. First is the combination of recent estimates 
of the Hubble parameter, which favour h = iifo/(100 km s~^ 
Mpc"^) > 0.55 (e.g., Tanvir et al. 1995; Bureau, Mould & 
Stavely-Smith 1996; Sandage et al. 1996; Riess, Press & Kir- 
shner 1996; Giovanelli 1996), with estimates of globular clus- 
ter that imply an age of the universe to ^ 12 Gyr (Chaboyer 
et al. 1996; Safaris, Degl'Innocenti & Weiss 1996; Renzini et 
al. 1996). These estimates are difficult to reconcile if f2o ~ 1- 
Second is the high baryon fraction in rich clusters of galaxies 
(e.g., Henry & Tucker 1979; Henriksen & Mushotzky 1985); 
this is incompatible with f2o = 1 if the ratio of baryons 
to dark matter in clusters is equal to the universal value 
and constraints on the baryon density derived from stan- 
dard nucleosynthesis theory and the observed light element 
abundances are correct (White et al. 1993b; Markevitch et 
al. 1996; Buote & Canizares 1996; Lubin et al. 1996). 

Traditional inflation models lead to a spatially flat uni- 
verse, and hence to f2o = 1 if the cosmological constant 
vanishes. The A-CDM models are perfectly compatible with 
standard inflation if one considers a non-zero cosmological 
constant to be physically reasonable. Open-bubble inflation 
models (Ratra & Peebles 1994,1995; Bucher, Goldhaber & 
Turok 1995; Yamamoto, Sasaki & Tanaka 1995) are vari- 
ants of the scenario suggested by Gott (1982) and Guth & 
Weinberg (1983), in which a spatially-open, inflating bubble 
nucleates via quantum tunneling within a spatially flat, ex- 
ponentially expanding (inflating), de Sitter spacetime. The 
post-nucleation inflation epoch stretches the open bubble to 
encompass the observable part of the universe, and the quan- 
tum fluctuations of the inflaton field inside this spatially- 
open bubble provide the seeds for structure formation. 

In the next section we describe our selection of models, 
with particular attention to the normalization of the linear 
power spectrum. We also describe our numerical simulation 
parameters. In Section 3 we compare the mass function of 
clusters in these simulations to observations. The models 
that appear viable after this comparison are the flat models 
with flo = 0.3 or 0.4, and the open models with flo — 0.4 
or 0.5. An alternative viable model is fio = 1 with a tilted 
(n < 1) primordial power spectrum. In Section 4 we present 
some statistical measures of the clustering of the nonlinear 
mass distributions in these models: the correlation function, 
the power spectrum, and the topology of isodensity surfaces. 

These analyses do not include any effects of biased or 
anti-biased galaxy formation, but they give a sense of what 
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sort of biases are needed to reconcile the model simulations 
with observations. In Section 5 we summarize our results 
and discuss the prospects for distinguishing these models 
with future observations. 



2 MODELS 

2.1 Parameter choices and normalization 

The primary parameter in the two classes of models we con- 
sider is the density parameter fio. However, the values of h 
and Q.B also have an influence on the shape of the matter 
power spectrum, and through this on the amplitude of mass 
fluctuations for the COBE-DMR, normalization. In light of 
recent evidence that tends to favour high values of h, we 
have decided to choose for each open model the value of h 
that gives to ~ 12 Gyr, i.e., the largest value that leaves the 
model marginally compatible with standard globular cluster 
age estimates. For each flat-A model we choose the value of 
h that gives to « 14 Gyr. For Oo = 1 models we take h = 0.5. 
Although our simulations do not include a separately treated 
baryon component, the value of the baryon fraction has a 
modest influence (weaker than the influence of h) on the 
shape and normalization of the power spectrum. To fix the 
initial amplitude and shape of the matter power spectrum, 
we adopt f^B ~ 0.0125/i~^, the value advocated by Walker et 
al. (1991) on the basis of the observed light element abun- 
dances and standard nucleosynthesis. Recent estimates of 
the deuterium abundance in high-redshift Lyman-limit ab- 
sorbers suggest that this value may be a factor of two high 
(Songaila et al. 1994; Carswell et al. 1994; Rugers & Hogan 
1996a,b; see however Tytler, Buries & Kirkman 1996) or 
low (Tytler, Fan & Buries 1996; Buries & Tytler 1996; see 
however Wampler 1996; Songaila, Wampler & Cowie 1996). 
On the scales modelled by the simulations in this paper, 
which are much smaller than present values of the Hubble 
or curvature radius, the shape of the energy-density per- 
turbation power spectrum is quite accurately determined 
by taking account of the usual effects of matter and radia- 
tion in processing the primordial energy- density perturba- 
tion power spectrum. We adopt the Bardeen et al. (1986) 
formula for the matter power spectrum, 



P{k) 



[1 + 3.89g -f (16.1q)2 + (5.46q)3 + (6.71^)'*]^''^ 
[ln(l + 2.34q)]' 



(2.34g)2 



where q — (fc/r)/iMpc ^. The shape parameter F is 



F = Q.i)hciqy(—Q,^ 



f^B/f^O 



(1) 



(2) 



(Sugiyama 1995). The slope of the primordial power spec- 
trum, n, we take to have the scale-invariant value (n = 1), 
except for the case of our one tilted model which has n ~ 0.8. 
This smooth analytic fit to the matter power spectrum is 



Note that lower values of tg might not be unreasonable. For 
example, Alcock et al. (1996) argue for to = 12 it 1.5 Gyr on the 
basis of a revised RR Lyrae distance scale. A smaller to allows for 
a larger h, hence more small-scale power and so lower Slo- 



fairly accurate in the limit r2B/fio ^ 1, but it does not re- 
produce the wiggles in the power spectrum caused by oscil- 
lations in the photon-baryon fluid when the baryon fraction 
is significant (Sugiyama 1995). We expect these wiggles to 
have minimal effect on cluster masses or small scale cluster- 
ing of the mass distribution, especially for the relatively high 
values of f2o that are favoured by our analysis in Section 3. 

The normalization of the matter power spectrum makes 
use of spectral information on the very large scales probed 
by the CO-B S-DMR experiment. For the f2o = 1 and open 
models, we take our values of ug, the rms, linearly extrapo- 
lated mass fluctuation in spheres of radius 8/i~^ Mpc, from 
Table 1 of Gorski et al. (1995), which is based on an analy- 
sis of the DMR two-year galactic-frame 53 and 90 GHz sky 
maps, with the observed quadrupole anisotropy moment ex- 
cluded from the analysis. For the flat models with a cosmo- 
logical constant, we take values of (Tg from Table 2 of Ra- 
tra et al. (1997), which is based on an analysis of the DMR 
two-year ecliptic-frame 31.5, 53, and 90 GHz sky maps, with 
the observed quadrupole anisotropy moment excluded from 
the analysis (Bunn & Sugiyama 1995). Similar values of erg 
for the flat-A models are obtained by Sugiyama (1995) and 
Stompor, Gorski & Banday (1995), and for the f2o = 1 
model by Bunn, Scott & White (1995), Sugiyama (1995), 
and Stompor et al. (1995). For the flat-A models these es- 
timates assume a scale-invariant primordial matter power 
spectrum (spectral index n = 1 in P(k) oc fc"), as predicted 
by the simplest inflation models. The simplest open-bubble 
inflation models, which we adopt for our open CDM mod- 
els in this paper, also produce matter power spectra that 
have this scale-invariant form on scales much smaller than 
the present Hubble length (Ratra & Peebles 1994). Except 
for the tilted fio = 1 model, our normalizations ignore the 
possible influence of primordial gravity waves or mild early 
reionization on the COB£'-DMR anisotropics. 

There are various ways to analyze the DMR two-year 
data that lead to slightly different erg normalizations for a 
specified cosmological model — one can, for example, use 
either the ecliptic- or galactic-frame sky maps and either in- 
clude or ignore the observed quadrupole anisotropy moment 
(which is the multipole most sensitive to the Galaxy model) 
in the analysis (e.g., Stompor et al. 1995). There is thus 
no unique central ug value for a "DMR-normalized" cosmo- 
logical model, but for our purposes the differences between 
these different methods are not very significant relative to 
the observational uncertainties in cluster masses discussed 
below in Section 3. 

After we had completed our simulations, results from 
the analyses of the COBE-T>MR, four-year data (Bennett 
et al. 1996; Gorski et al. 1996a, and references therein) be- 
came available. In addition to new data, these analyses in- 
corporate more sophisticated approaches to removing Milky 
Way contamination. The net effect is a slight downward 
shift in central estimates of the power spectrum normal- 
ization for the models, mostly a consequence of the new 
(four-year DMR) Milky Way exclusion cut (Banday et al. 
1997). Different combinations of anisotropy maps, treatment 
of the quadrupole, and corrections for Milky Way emission 
again lead to slightly different estimates of erg (Gorski et 
al. 1996a,b), but relative to the average of these DMR four- 
year estimates our adopted normalizations are typically high 
by 6 — 12%, on the same order as the l-cr DMR four-year 
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Table 1. Simulation Parameters 



Model 


Qo 


Ao 


h 


to/Gyr 


f^B 


r 


0"8 


a} 


Afstops 


01 


0.1 


0.0 


0.75 


11.7 


0.022 


0.05 


0.1 


0.98 


2 


02 


0.2 


0.0 


0.70 


11.8 


0.026 


0.12 


0.3 


0.28 


35 


03 


0.3 


0.0 


0.65 


12.2 


0.030 


0.17 


0.5 


0.15 


93 


04 


0.4 


0.0 


0.65 


11.7 


0.030 


0.23 


0.75 


0.1 


168 


05 


0.5 


0.0 


0.6 


12.3 


0.035 


0.27 


0.9 


0.08 


254 


LI 


0.1 


0.9 


0.9 


13.9 


0.035 


0.07 


0.7 


0.15 


150 


L2 


0.2 


0.8 


0.75 


14.0 


0.022 


0.13 


0.9 


0.12 


220 


L3 


0.3 


0.7 


0.65 


14.5 


0.030 


0.17 


1.05 


0.101 


266 


L4 


0.4 


0.6 


0.6 


14.5 


0.035 


0.21 


1.1 


0.09 


275 


L5 


0.5 


0.5 


0.6 


13.5 


0.035 


0.27 


1.3 


0.07 


331 


El 


1.0 


0.0 


0.5 


13.1 


— 


0.5 


1.35 


0.0605 


327 


E2 


1.0 


0.0 


0.5 


13.1 


0.05 


— 


0.55 


0.2 


200 



data uncertainty. Depending on spectrum and parameter 
values, the DMR four-year data uncertainty (including all 
known systematic uncertainties) in the ag normalization is 
~ ±(16 - 19)% (2-cr) (Gorski et al. 1996b). 

In addition to this series of models based on the simplest 
inflation power spectra, we consider an f2o = 1 model with a 
tilted primordial power spectrum, n < 1. We have chosen the 
degree of tilt so that when the (power-law inflation model) 
gravity wave contribution to the DMR CMB fluctuations is 
included (Bunn & White 1997; Liddle et al. 1996), the DMR- 
normalized ag = 0.55, which is about the value that White et 
al. (1993a) conclude is needed to explain the mass function 
of galaxy clusters in an f2o = 1 universe. This prescription 
and a Hubble parameter of ft = 0.5 leads to n = 0.803, for 
which the reduction in the implied ag, caused by the inclusion 
of gravitational waves is 27%. 

The values of f2o, Ao (in units of 3Hq), h, to, Q,b, T, 
and as for all of our models are listed in Table 1. We label 
the five flat-A models L1-L5, the five open models 01-05, 
the DMR- normalized Einstein-de Sitter (fio = 1) model El, 
and the tilted Einstein-de Sitter model E2. Note that for the 
model El r = 0.5, as is conventional in the standard CDM 
model, so here we have effectively assumed fls = 0. 



2.2 Simulations 

Each of our simulations models a periodic cube of side 
345. 6/i^^ Mpc (comoving). We want large boxes so that 
we can get accurate statistical estimates of cluster abun- 
dances, and also because, in future work, we will use these 
simulations to create mock redshift catalogues of the Sloan 
Digital Sky Survey (see Gunn & Weinberg 1995) and the 
Anglo- Australian 2-Degree Field Survey (see CoUess 1996). 
We evolve the simulations using a modified version of 
Hugh Couchman's Adaptive Particle-Particle-Particle-Mesh 
(AP^M, Couchman 1991) N-body code. The P^M algorithm 
(Efstathiou et al. 1985) allows us to simulate large volumes 
with a force and mass resolution sufficient to yield accurate 
cluster masses. 

Instead of starting our simulations from a regular 
particle grid, we use the technique described by Baugh, 
Gaztafiaga & Efstathiou (1995) and White (1996) to gen- 
erate "glass" initial conditions. We lay down a Poisson dis- 
tribution of 192'' particles, then integrate these with the N- 



body code and sign-reversed gravitational forces until they 
reach a state at which the gravitational forces on all parti- 
cles practically vanish. With this approach, the initial par- 
ticle distribution is not regular, but random fluctuations in 
the particle density field do not seed the growth of spuri- 
ous structures. For each simulation, we create a Gaussian 
random density field on a 192^ grid, using the same Fourier 
phases from one model to the next but varying the power 
spectrum in accord with the model parameters. We apply 
the Zel'dovich approximation to this density field to com- 
pute displacements and peculiar velocities for the particles, 
interpolating from the grid to the particle positions. Simu- 
lations with glass initial conditions and grid initial condi- 
tions yield similar statistical results once they are evolved 
into the nonlinear regime (Baugh et al. 1995; White 1996), 
but simulations with glass initial conditions do not retain a 
grid signature in uncoUapsed regions. We set the softening 
parameter of AP^M's triangular-shaped cloud force law to 
77 — 270h^^ kpc, 0.15 of the grid spacing. This choice cor- 
responds approximately to a gravitational softening length 
e = ri/3 — 90h^^ kpc for a Plummer force law, and we adopt 
e as our nominal force resolution. 

The initial expansion factors of the simulations a' , listed 
in Table 1, were determined by setting the amplitude of the 
initial power spectrum at the Nyquist frequency of the par- 
ticle grid to be 0.3^ times that for an equivalent Poisson 
distribution of particles. Thus Pinitiai(fcN) ~ 0.3^/n, where 
n is the mean particle density and the Nyquist frequency is 
/cn = 7rn^" — (27r/3.6)/i Mpc~^. The residual power in the 
glass configuration is only 0.5% of that in a Poisson distri- 
bution at the Nyquist frequency and drops very rapidly at 
longer wavelengths (see figure A2 of Baugh et al. 1995). Thus 
this choice is safely in the regime where the initial density 
fluctuations are large compared to those present in the glass 
and yet in the linear regime where the Zel'dovich approx- 
imation remains accurate. In particular, no shell-crossing 
has occurred. (The Ol model has an amplitude lower than 
this choice at the final time. Consequently, we started this 
simulation with fluctuation amplitude a factor of 2 lower.) 
The size of the time step Aa was chosen so that the fol- 
lowing two constraints were satisfied throughout the evolu- 
tion of the particle distribution. First, the rms displacement 
of particles in one timestep was less than rj/A. Second, the 
fastest moving particle moved less than rj in one timestep. 
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Figure 1. The evolved particle distribution in a 20h^^ Mpc thick slice of the Qq = 0.4, open universe simulation. The large panel shows 
the full cross section of the simulation box, while the two lower panels show successively expanded views, lOOh^^ Mpc and 20h^^ Mpc 
on a side respectively. 
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LI: Oo=0.1 Aa=0.9 
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Figure 2. Evolved particle distributions from the six DMR-normalized, flat universe models. The slices are 20h Mpc thick and 
100/i~^ Mpc on a side, and they show the same portion of the simulation volume that is shown in the lower left panel of Fig. 1. 



Initially these two constraints are comparable, but at late 
times the latter constraint is more stringent, particularly in 
low Q.0 simulations. We monitored energy conservation us- 
ing the Layzer-Irvine equation (equation 12b of Efstathiou 
et al. 1985) and found that for this choice of timestep en- 
ergy conservation with a fractional accuracy of better than 
0.3% was achieved. We also tested the inaccuracy incurred 
by these choices of starting amplitude and timestep by com- 



paring the final particle positions with two additional fio = 1 
simulations that were run starting from a fluctuation ampli- 
tude a factor of two lower and using time steps a factor of 
two smaller. In each case we found the final particle posi- 
tions to agree very accurately, with rms differences of less 
than e. More importantly, the correlation functions of each 
particle distribution were indistinguishable at scales larger 
than e — 90/i^^ kpc. Thus the statistical clustering proper- 
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Figure 3. Same as Fig. 2 but for open universe models. The bottom right panel shows the tilted n = 0.803, f2o = 1 model, E2. 



ties of these simulations have their resolution limited by the 
particle mass and force softening and not by the choice of 
timestep or starting redshift. 

The large panel of Fig. 1 shows a slice 20h~^ Mpc thick 
(6% of the simulation volume) through the evolved parti- 
cle distribution of the flo = 0.4, open universe model (04). 
The distribution follows the familiar pattern of elongated, 
sharp, high density features and rounded voids that arises 
generically from the action of gravitational instability on 



Gaussian initial conditions. To illustrate the dynamic range 
of the simulation, we show expanded regions lOO/i^^ Mpc 
and 20h~^ Mpc on a side in the lower two panels. At the 
resolution level of the 20h~^ Mpc panel, one sees that the 
coherent large scale features are composed mainly of dense, 
ellipsoidal clumps, connected by a sprinkling of isolated par- 
ticles. 

Fig. 2 shows lOO/i"^ Mpc x lOOft"^ Mpc x 20h'^ Mpc 
slices from our six DMR-normalized, fiat universe models, 
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L1-L5 and El. Moving from f2o = 0.1 through successively 
higher values of fio, we see a steadily increasing degree of 
clustering, primarily because the value of o» implied by 
DMR- normalization is an increasing function of f2o- The 
mass per particle also grows in proportion to fio, so the 
masses of the collapsed structures increase more rapidly 
with fio- The ag, normalization of the f2o = 1 model is 
only slightly larger than that of the f2o = 0.5 model, but 
its structure is choppier because its initial conditions have 
more small scale power (F = 0.5 vs. T = 0.27). 

Fig. 3 shows slices from the five open models, 01-05, 
and from the tilted fio = 1 model, E2. We see the same 
trends as in Fig. 2, but here they are much stronger, because 
in open models the erg normalization implied by the DMR 
data is more sensitive to the value of f2o- The f2o = 0.1 
model has virtually no collapsed structure on scales that 
our simulation can resolve, and structure in the fio = 0.2 
and fio = 0.3 models is rather anemic. The structure in the 
E2 model is intermediate between that of the 03 and 04 
models, which have bracketing values of ug- 

From these plots we can already anticipate the qual- 
itative results for cluster masses. In each class of models, 
cluster masses will increase with increasing f2o because of 
the greater level of clustering and the higher mass density. 
The dependence in the open models will be stronger than in 
the flat models. 



3 THE CLUSTER MASS FUNCTION 

The space density of massive clusters is a sensitive function 
of fio and os ■ Since these two parameters vary in concert for 
DMR-normalized CDM models, the cluster abundance can 
provide a tight constraint on fio- More important still, the 
abundance of clusters as a function of mass can be mea- 
sured directly from N-body simulations in a way that is 
insensitive to uncertain assumptions about galaxy forma- 
tion. We must assume that galaxies form efficiently enough 
in cluster-mass potential wells that these are indeed identi- 
fied as galaxy clusters, but this seems more plausible than 
the alternative assumption that the universe contains many 
cluster-mass dark halos that harbour few galaxies and lit- 
tle X-ray emitting gas. Most of the uncertainty in using the 
cluster abundance as a cosmological test comes from the sys- 
tematic uncertainties in observational estimates of cluster 
masses, a point to which we will return later in this Section. 
The predicted present-day abundance of clusters in our 
DMR-normalized models is given in Fig. 4. The top panel 
shows the abundance of clusters as a function of the mass 
contained within the Abell radius, 7?a = 1.5/i^^ Mpc, of 
each cluster. Following White et al. (1993a) we also plot, 
in the lower panel, the abundance as function of the mass 
within -Ra/3, since this can often be more robustly esti- 
mated from observations (e.g., van Haarlem, Frenk & White 
1997). In both cases cluster centres were identified using a 
friends-of-friends group finder with a small value of the link- 
ing length, 0.17/n^'^ = 306/i^^ kpc. This value picks out the 
high density knots of the mass distribution with overdensity 
of approximately 1000. Masses were then measured within 
spheres centered on the centre-of-mass of these knots. Dou- 
ble counting was eliminated by ranking the cluster centres 
by the mass of the friends-of-friends group and assigning 



particles to the most massive within Ra- This is a small 
correction, and the results are very similar if one simply ex- 
cludes clusters that lie within _Ra of a more massive cluster. 

In all models, the abundance of clusters declines rapidly 
with increasing mass. Clusters are, on average, more massive 
in models with larger Q.q. The number density of clusters of a 
given mass increases rapidly with Oq (especially for low Oo), 
and many more clusters are formed in fiat universes than in 
open universes with the same value of r2o. For example, there 
arc ~ 6 times fewer clusters with M{Ra) > 3 x IQ^'^h'^ M© 
in the flat ilo = 0.4 model and ~ 25 times fewer in the open 
0() = 0.4 model than there are in the Oo = 1 El model. 
Flat models with f2o < 0.1 and open models with fio < 0.3 
produce a negligible abundance of rich clusters. 

The cluster abundance has been applied as a cosmo- 
logical test by, among others. White et al. (1993a), Viana 
& Liddle (1996), Eke et al. (1996b), Oukbir & Blanchard 
(1997), Oukbir, Bartlett, & Blanchard (1997), Kitayama & 
Suto (1996), and Gorski et al. (1996b). These authors com- 
pared model predictions based on the Press & Schechter 
(1974) formalism (or generalizations thereof) with the ob- 
served abundance of clusters measured as a function of mass 
in the first study and as a function of X-ray temperature in 
the later studies. T The advantage of using a grid of N-body 
simulations over the more fiexible analytic approach is that 
cluster masses can be directly measured in the simulations, 
thus bypassing the need to assume a model relating pre- 
dicted cluster masses to either measured masses or X-ray 
temperatures. 

To compare the predictions in Fig. 4 to observations re- 
quires a measurement of the abundance of clusters in a spec- 
ified mass range. The open circles in the figure show Bahcall 
and Cen's (1993) analytic fit to their estimated mass func- 
tion of groups and clusters in the local Universe. In agree- 
ment with their conclusions, we find that a fiat model with 
Q.{) ~ 0.3 provides a good match to their data, whereas a 
DMR-normalized Einstein-de-Sitter model with fiducial pa- 
rameter values (Model El) overpredicts the abundance of 
clusters by a large factor. These conclusions, however, are 
uncertain because of the difficulties inherent in defining a 
complete sample of clusters over a large range of masses. 
Further systematic uncertainties could be introduced by the 
heterogenous mass estimators employed by Bahcall & Cen. 
[We note, however, that, given the uncertainties, Bahcall & 
Cen's mass function is in fairly good agreement with the 
temperature function derived by Eke et al. (1996b) from 
Henry & Arnaud's (1991) X-ray data and a simple model 
relating cluster masses and X-ray temperatures.] A more 
reliable comparison between models and observational data 
may be made by considering a single robust statistic charac- 
terising the cluster abundance, rather than the mass func- 
tion as a whole. The median mass of Abell clusters with 



T The first four sets of authors expressed the requirement that a 
model should reproduce the observed cluster abundance in terms 
of an acceptable range of values for erg as a function of Qq ■ Com- 
parisons of these with the values of erg required by the DMR 
normalization are given, for example, in figure 13 of Eke et al. 
(1996b), figure 4 of Viana & Liddle (1996) and tables 9-12 of 
Gorski et al. (1996b). On the other hand, Kitayama & Suto (1996) 
computed cluster temperature and luminosity functions for DMR- 
normalized models and compared these to the observations. 
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Figure 4. Cumulative cluster mass functions; n{M) is the number density of clusters whose mass exceeds M . In the upper panel, masses 
are defined within the Abell radius Ra = 1.5/i~^ Mpc. In the lower panel, masses are defined within Ra/3. In each panel, the upper 
dashed line represents the DMR- normalized fiducial Uq = 1 model. El, and the lower dashed line the DMR-normalizcd tilted Qq = 1 
model, E2. Dotted lines represent the A-dominated models L1-L5, with Qq increasing from bottom to top. Solid lines represent the open 
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richness class R> 1, advocated by White et al. (1993a), has 
been shown to be robust to the distortion of the shape of 
the measured cluster mass function caused by line-of-sight 
projections (van Haarlem et al. 1997). The abundance of 
R>1 clusters, Na = 8 x 10~''/i^Mpc~^, is weU established 
from X-ray and optical studies (Bahcall & Soneira 1983; 
Scaramella et al. 1991; Efstathiou et al. 1992b). The results 
of the comparison do not depend sensitively on the adopted 
value of Na because the predicted number density of clus- 
ters declines very steeply with mass. For the same reason, 
however, the results are quite sensitive to the adopted value 
of the median mass. 

White et al. (1993a) give the range McIus{Ra) = (4.2 - 
5.5) X lO^'^/i^^MQ for the median mass within the Abell ra- 



dius of i? > 1 Abell clusters. The lower limit comes from es- 
timates using X-ray data and the upper limit from estimates 
using velocity dispersions of cluster galaxies. As White et al. 
emphasize, the statistical significance of this range is diffi- 
cult to quantify because most of the uncertainties arise from 
systematic errors in the mass determinations. The quoted 
range is meant to take at least some of these systematics 
into account. 

We adopt a slight modification of this range. We retain 
the upper mass limit advocated by White et al. (1993a), but 
reduce the lower mass limit to 3.0 x lO^'^/i^^M©. This lower 
value follows from two considerations. First, the reanalysis 
by Eke et al. (1996b) of the X-ray data compiled by Henry 
& Arnaud (1991) yields an X-ray temperature 3.3 keV for 
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a median R = 1 cluster rather than the 3.6 keV adopted 
by White et aL (1993a). Second, White et al. extrapolated 
the cluster mass from the virial radius to Ra assuming an 
r^^ density profile. This extrapolation amounts to a factor 
of 1.35 in mass. If cluster density profiles are significantly 
steeper at the virial radius, as suggested by Navarro, Frenk 
& White (1995), then a more modest factor is appropriate. 
Hence we adopt 3.0 x lO^''/i~^M0 as a conservative lower 
limit. For masses within Ra/3, we simply reduce the mass 
range by a factor of three at the upper and lower ends, which 
is correct if cluster density profiles are oc r~^. In practice, 
the estimated masses within Ra/3 are probably more robust 
than those within Ra, since the X-ray data typically do 
not extend all the way to 1.5h~^ Mpc. However, we reach 
similar conclusions about the viability of models using either 
measure. 

In addition to the traditional methods for estimating 
cluster masses, based on the kinematics of their galaxies 
and/or on the hydrodynamics of their X-ray emitting gas, 
a relatively new method, based on the gravitational lensing 
properties of the cluster dark matter, has recently been ap- 
plied to several clusters. The extent to which all three meth- 
ods give consistent results is still a matter of debate. All of 
them require a number of assumptions of varying degrees of 
plausibility (see, e.g.. White et al. 1993a). With standard as- 
sumptions, estimates based on velocity dispersions often give 
somewhat larger masses than estimates based on X-ray data. 
There are examples for which the mass inferred from grav- 
itational lensing differs significantly from the mass inferred 
from X-ray or dynamical data (e.g., Fahlman et al. 1994) 
and others where the agreement is good (e.g.. Squires et al. 
1996). At present, the sample of clusters analyzed by grav- 
itational lensing is small, and it consists predominantly of 
clusters at intermediate redshift, where the traditional meth- 
ods are particularly uncertain. An important argument for 
the robustness of dynamical estimates based on galaxy mo- 
tions is the fact that blue and red galaxy populations yield 
compatible masses, despite having different density and ve- 
locity dispersion profiles (Carlberg et al. 1997c). 

From Fig. 4, we see that the DMR-normalized CDM 
models that produce an acceptable cluster abundance are 
the open 04 and 05 models, the flat L3 and L4 models, 
and, by construction, the tilted flo — 1 model E2. Thus we 
estimate that fio > 0.25 is required to produce the observed 
abundance of rich clusters in the spatially-flat models and 
that this lower limit rises to f2o > 0.4 for open models. The 
corresponding upper limits are flo < 0.4 for the flat models 
and Qo ^ 0.5 for open models. The fio = 1 model with 
the scale-invariant n = 1 primordial power spectrum. El, 
produces an abundance of 7? > 1 clusters that is 20 times 
too high. 

These constraints are in accord with more approximate 
ones derived earlier (e.g.. Eke et al. 1996b; Kitayama & Suto 
1996; Gorski et al. 1996b). Note, however, that our analysis 
does not account for the known uncertainties in estimates of 
to, ^B, and the DMR normalization. These uncertainties will 
somewhat broaden the above allowable ranges of Qq. For ex- 
ample, Gorski et al. (1996b) compare DMR-normalized val- 
ues of (78 to those deemed necessary by Eke et al. (1996b) to 
explain the observed cluster temperature function, and ac- 
counting for the 2-(j uncertainty in the DMR normalization 
and some of the uncertainty in to and f^B, they find an al- 



lowable range of 0.3 < Qq < 0.6 (~2-ct) 
infiation model. 



for the open-bubble 



4 OTHER MEASURES OF CLUSTERING 

We now compute some statistical characteristics of the real 
space mass distribution for the five models from our suite 
that produce observationally acceptable cluster abundances. 
The statistics of the galaxy distributions in these models will 
depend on the relation between galaxies and mass ("bias" 
or "antibias"), and, if they are measured in redshift space, 
on the distortions caused by peculiar motions. The measures 
presented here will give some sense of how the models com- 
pare to each other and indicate what kinds of bias would 
be needed to obtain acceptable fits of our N-body results to 
observational data. 

Figs. 5 and 6 show the mass autocorrelation functions 
^(r) and power spectra P{k) from these five models. We 
also show estimates of the real space galaxy autocorrela- 
tion function and power spectrum derived by Baugh (1996) 
and Baugh & Efstathiou (1993) from the angular correlation 
function of the APM galaxy catalogue (Maddox et al. 1990). 
Since the APM catalogue is fairly deep, the inferred (,{r) and 
P{k) depend somewhat on the assumed redshift evolution of 
clustering. We show Baugh's and Baugh and Efstathiou's re- 
sults for two different assumptions, linear theory evolution 
for fio = 1 (open squares) and clustering fixed in comov- 
ing coordinates (filled squares). These assumptions probably 
bracket the true evolution, and the difference between them 
illustrates the systematic uncertainty associated with this 
effect. The values crg.gai, the rms fluctuation of the density 
of galaxies in spheres of radius 8h~^ Mpc, implied by these 
two results are (Tg^gai ~ 0.82 (fixed in comoving coordinates) 
and (78, gal ~ 0.92 (linear evolution). The inferred clustering 
amplitude, which we can characterize by the value of (Js.gai, 
is also subject to other systematic uncertainties. Uncertain- 
ties in the APM selection function and in the assumed cos- 
mological model both introduce uncertainties at the level 
of a few percent. Furthermore, systematic errors exist in 
APM angular correlation function at some level. Maddox, 
Efstathiou & Sutherland (1996) analyse a wide range of pos- 
sible systematic errors and apply corrections for recognized 
problems with star-galaxy classification at faint magnitudes 
and other small systematic errors. These corrections lead to 
o"8,gai ~ 0.96 for the same cosmological model and redshift 
distribution assumed above with a clustering evolution inter- 
mediate to the two cases above. The error bars shown on the 
APM data in Figs. 5 and 6 indicate the statistical errors and 
are estimated from the dispersion among four nearly equal 
zones of the APM survey. For the most part, the sources of 
systematic error discussed above tend to alter the amplitude 
of ^(r) and P{k) without affecting their shape. However, the 
points at largest separation (r > 40/i^^ Mpc) and lowest 
wavenumber {k < 0.06ftMpc^^) probably have uncertain- 
ties that are dominated by systematic errors (Maddox et al. 
1996). 

Fig. 5 shows the mass autocorrelation functions ^(r) 
from the five acceptable models, measured from twice the 
effective gravitational softening length, 2e = 180h~^ kpc, 
out to 35/i~^ Mpc, ~ 1/10 of the simulation box size. Fig. 6 
shows the mass power spectra of the five models. On the 
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Figure 5. Mass autocorrelation functions for the five models from our suite of N-body simulations that yield acceptable cluster masses. 
Short-dashed and dotted lines show the fiat models with f!o = 0.4 and 0.3, respectively, long-dashed and dot-dashed lines show the open 
models with Q,q = 0.5 and 0.4, respectively, and the solid line shows the tilted Qq = 1 model. Results are shown from r = 35/i~^ Mpc, 
1/10 of the box size, down to r = 0.2/i~^ Mpc, about twice the effective gravitational softening length e = 0.09/i~^ Mpc. Solid squares 
with error bars show the real space galaxy correlation function inferred by Baugh (1996) from the APM angular correlation data assuming 
that it remains fixed in comoving coordinates. Open squares show the galaxy correlation function derived assuming linear theory evolution 
of the correlation function with redshift — statistical errors on these points (not shown) are similar to those on the solid points. 



largest scales the linear theory power spectra are plotted 
and these are then continued to higher k using estimates 
from the N-body simulations. We measure the model power 
spectra in the N-body distributions by FFT after cloud-in- 
cell (CIC) assignment of the particle distribution onto a 192^ 
grid. We correct the FFT power spectrum for the effect of 
CIC convolution. In two cases (models E2 and L3), the linear 
theory power spectrum is shown for all k so that the effect 
of non-linear evolution on P{k) can be seen. Error bars are 
shown for model E2, which illustrate the uncertainty in esti- 
mating P{k) on very large scales where there are few modes 
in the simulated cube. The relative biases required in the 
models if they are to reproduce the observed galaxy cluster- 
ing are shown most clearly in Fig. 7, which shows the ratio 
of APM galaxy power spectrum to the model mass power 
spectra. 

The L4 and L3 models, which have (F, as) = (0.21, 1.1) 
and (0.17, 1.05), respectively, have mass correlation func- 



tions which trace the APM galaxy correlations fairly well at 
large separations, r > lOfe"^ Mpc. However, as previously 
noted (e.g., Efstathiou et al. 1990; Klypin et al. 1996; Mad- 
dox et al. 1996; Peacock 1996), their mass correlation func- 
tions steepen considerably at r < 4fe^^ Mpc, rising above 
the APM galaxy data by a factor --^ 3 at r = \h^^ Mpc. To 
match the observed galaxy correlation function, these mod- 
els would require a strong antibias between galaxies and 
mass on these scales but could not tolerate much antibias 
or bias at large scales. The 05 model has a higher F (0.27) 
and a lower ag (0.9). Its mass correlation function traces the 
APM galaxy data for scales where ^(r) ~ 1. The mass cor- 
relation function again steepens at smaller scales and rises 
above the APM galaxy data, though not as much as for the 
A models. It falls below the APM data at scales larger than 
r ~ lQh~^ Mpc. To fit the data, this model would require a 
positive bias at large scales and an anti-bias at small scales. 
The 04 model, with (F, erg) = (0.23,0.75), traces the APM 
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Figure 6. Mass power spectra for the five models shown in Fig. 5, with the same correspondence between model and line type. For each 
model the heavy curve shows the results of the N-body simulation and the light curve the corresponding linear theory power spectrum. 
For the model E2, we indicate the \-a errors due to the finite number of modes in the simulation cube. The smallest fc value in the E2 
N-body P(}i) is the fundamental frequency of the 345.6/i~^ Mpc simulation box; the largest is O.Sfejv, where fcjv is the Nyquist frequency 
of our 192^ FFT grid. For the other models we omit the error bars and do not plot the first three, noisy data points. For models E2 
and L3 the linear theory curves are plotted over the full range of k to illustrate the effect of non-linearities at small scales. For the other 
models the linear theory power spectra are just plotted as an accurate extrapolation of the model power spectra to small fc. The filled 
squares with error bars show the real space galaxy power spectrum inferred by Baugh & Efstathiou (1993) from the APM angular power 
spectrum, assuming that clustering remains fixed in comoving coordinates. Open squares (without error bars) show the corresponding 
galaxy power spectrum assuming linear theory evolution of clustering with redshift for an S7() = 1 universe. 



galaxy data fairly well for r < 5/i^^ Mpc. Because mass 
fluctuations are less nonlinear in this model, the mass corre- 
lation function does not show the marked steepening inside 
the correlation length seen in the earlier models. However, 
the correlation function falls below the APM data at large 
separations, so the 04 model would require positive bias on 
these scales in order to match the data. The E2 model, with 
(r, erg) = (0.45, 0.55) and n — 0.803, requires a positive bias 
on all scales to match the APM galaxy data, as expected. 
Again the low amplitude means that there is no marked 
steepening of the slope inside the correlation length. The 
required bias is roughly independent of scale except at the 
largest separations, where the E2 correlation function curls 
away from the data. The APM data are well fit on these 
scales by a model with an n = 1 primordial power spectrum 
and r ~ 0.15 - 0.2. The E2 model fits better than an n = 1 



model with the same value of F = 0.45, but it nonetheless 
yields smaller correlations than the APM data on scales of 
30 - lOOh-^ Mpc. 

Overall, the power spectrum comparison confirms our 
conclusions from the ^(r) comparison (as expected, since 
the two measures are a Fourier transform pair). The E2 
model requires a positive bias on all scales to match the 
APM galaxy data; the 04 and 05 models require a positive 
bias on large scales (fc ~ 0.05 h Mpc^^), with 05 also re- 
quiring antibias on small scales; and the L3 and L4 models 
require a weak antibias on large scales (previously noted in 
linear theory mass P{k) and CMB anisotropy computations, 
e.g., Stompor et al. 1995; Ostriker & Steinhardt 1995; Ra- 
tra et al. 1997; Ganga et al. 1996c) and a strong antibias on 
small scales. 

The correlation function and the power spectrum quan- 
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Figure 7. The ratio of the APM galaxy power spectrum, as inferred by Baugh & Efstathiou (1993) for clustering fixed in comoving 
coordinates, to the mass power spectrum of each of the five models shown in Fig. 5 and 6. The correspondence between model and line 
type is the same as in the previous two figures. The error bars (for clarity only plotted on two of the models) indicate only the statistical 
error arising from the APM galaxy power spectrum. Uncertainties associated with the assumed evolution of clustering and APM selection 
function (see text) can lead to a systematic shift of ^ 0.2 in \og{PAPM(k)/P{k)). Also, for k < 0.06/iMpc~ , systematic errors in the 
APM angular correlation function probably result in errors in log(P^pM(fe)/f (fc)) which arc larger than the plotted statistical errors 
(Maddox et al. 1996). 



tify the rms amplitude of fluctuations as a function of scale, 
but they contain no further information about the global 
structure of the density field. As a complementary statis- 
tic, we consider the topology of isodensity contour sur- 
faces, following the methodology of Gott, Melott, & Dick- 
inson (1986) and Gott, Weinberg, & Melott (1987). We 
smooth the mass density field by convolution with a Gaus- 
sian filter of radius Ts = 4.5/i~^ Mpc (corresponding to 
A = 4.5^ = 6.4/i"^ Mpc with the smoothing filter defi- 
nition of Gott et al. 1987), then measure the genus Gs of 
isodensity contours at a range of threshold densities. Fig. 8 
plots the genus per unit volume, g = G's/(345.6/i~ Mpc) , 
multiplied by the scaling factor An'^r'i., as a function of v, 
which is defined implicitly in terms of the contour's enclosed 
fractional volume / through the equation 



/ = 



2/2 



dx. 



(3) 



An advantage of the volume-weighted procedure implied by 
equation (3) is that it makes the "genus curve" insensitive 
to biasing or antibiasing of the galaxy distribution, in con- 
trast to ^(r) or P{k). To the extent that biasing preserves a 
monotonic relation between the smoothed mass density and 
smoothed galaxy density, it does not alter the genus curve, 
even if the relation is nonlinear. In practice, plausible bias- 
ing models can slightly distort the genus curve (Park & Gott 
1991), e.g., by attenuating the walls of galaxies that separate 
underdense regions, but we expect that the genus curves of 
galaxy density fields in these models would be close to the 
genus curves of the mass density fields shown in Fig. 8. 

A Gaussian random field with power spectrum P{k) ex 
fc" has a mean genus curve 

3/2 



3 + n 
6 



(l-^^e- 



V2 



(4) 



(Doroshkevich 1970; Adler 1981; Bardeen et al. 1986; Hamil- 
ton et al. 1986). With the volume- weighted convention for 
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Figure 8. Topology of isodcnsity contour surfaces for the five models shown in Fig. 5. The topology is measured from mass density 



fields that have been smoothed with a Gaussian filter of radius Ts 
against the value of v defined in eq. (3). 



: 4.5/i Mpc. The scaled genus-per-unit-volume, 47r r^g, is plotted 



labelling contours, mildly nonlinear gravitational evolution 
has only a small eflect on the shape of the genus curve, 
though it does tend to lower the overall amplitude (Melott 
et al. 1988; Park & Gott 1991). The symmetric, W-shaped 
form of the curves in Fig. 8 reflects the Gaussian nature 
of the initial fluctuations in all of these models. The dif- 
ferent amplitudes reflect the differences in the logarithmic 
slope of the primordial power spectrum on the 4.5/i^^ Mpc 
smoothing scale (eq. [4]) and the differing degrees of nonlin- 
ear dynamical evolution corresponding to different values of 
an. 

A disadvantage of the topology statistic is that mea- 
surements from "small" redshift samples (e.g., flux-limited 
redshift surveys of a few thousand galaxies) suffer significant 
systematic biases because of the finite survey volume and 
the unknown boundary conditions for defining the smoothed 
density field (see, e.g., the discussion by Protogeros & Wein- 
berg 1997). These can be taken into account by comparing 
observational results directly to those of artificial catalogues 
drawn from N-body simulations, but this level of detailed 
modelling is beyond the scope of this paper. We have not 
included observational data in Fig. 8 because there is no sim- 
ple way to incorporate the effect of these sample-dependent 
systematic biases, but studies comparing mock catalogues to 



a number of galaxy redshift surveys show that the observed 
topology is generally consistent with models that have Gaus- 
sian initial fiuctuations with a F « 0.25 power spectrum 
(Gott et al. 1989; Park et al. 1992; Moore et al. 1992; Vo- 
geley et al. 1994; Protogeros & Weinberg 1997). Current 
topology data are probably not sufficient to distinguish be- 
tween the models illustrated in Fig. 8, but future redshift 
surveys like the 2dF and SDSS (CoUess 1996; Gunn & Wein- 
berg 1995) should yield precise topology measurements with 
minimal systematic errors that can easily distinguish most 
of these models from each other (the exception being 04 and 
05, which have nearly identical genus curves). Comparisons 
with high-precision data will also need to incorporate the 
impact of galaxy biasing, but we expect this to be relatively 
small for the reasons discussed earlier. 



5 DISCUSSION 

We have used large N-body simulations of CDM models 
to constrain the range of model parameters which, given 
some simple assumptions, are consistent with the cosmic mi- 
crowave background anisotropies measured by CO-B S-DMR 
and the abundance of rich galaxy clusters in the local uni- 
verse. We have then investigated in detail some basic cluster- 
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ing statistics for the dark matter in a selection of acceptable 
models. 

The combination of DMR anisotropy and cluster abun- 
dance provides an appealing way to constrain the parame- 
ters of a CDM universe because, by avoiding explicit refer- 
ence to galaxies, it bypasses the highly uncertain connection 
between the large-scale distributions of galaxies and mass. 
Nevertheless, this test does have a number of important 
limitations. The extrapolation of the fluctuation amplitude 
from the DMR scale to the cluster scale, r ~ 8/i^^ Mpc, 
requires an assumption about the shape of the primordial 
mass power spectrum emerging from inflation, e.g., a scale- 
invariant (n = 1) spectrum in flat models or the simplest 
open-bubble inflation spectrum (which is close to n = 1 on 
scales much smaller than the present Hubble length). This 
shape is not tightly constrained by the DMR data, and be- 
cause of the large lever arm between the DMR scale and 
8/i~^ Mpc, small changes in the spectral shape can produce 
large changes in the inferred value of ag. This uncertainty 
is compounded by the fact that in standard flat inflation 
models even small deviations from n = 1 can result in a 
significant contamination of the measured anisotropy signal 
by tensor modes, which produce gravity waves rather than 
density fluctuations (e.g., Salopek 1992; Crittenden et al. 
1993). 

The determination of erg from the observed cluster 
abundance has the great advantage of being essentially in- 
dependent of the shape of the power spectrum iiflo — 1 and 
of depending only very weakly on it if f2o < 1. The depen- 
dence on Ao is also weak; for flo > 0.2, the difference in as 
between flat and open models with the same value of Qq is 
always less than 10 per cent (Eke et al. 1996b). The main 
uncertainty stems from the difficulty in estimating cluster 
masses reliably. The expected number of clusters declines 
very steeply with cluster mass and, as a result, the estimate 
of (78, while not requiring a particularly accurate measure- 
ment of the abundance of clusters, is very sensitive to the 
estimates of cluster mass. In this paper, we have adopted 
a slightly wider range of cluster masses than than quoted 
by White et al. (1993a) for the median mass of i? > 1 
Abell clusters within the Abell radius Ra = 1.5h~^ Mpc: 
-McIus(-Ra) = (3.0- 5.5)/i"^M0. The lower end of this range 
corresponds to our revised estimate from X-ray data and the 
upper end to the estimate from dynamical data. 

In our series of low-f2o CDM simulations, we assume the 
simplest primordial inflation power spectra, a Hubble con- 
stant that yields to « 14 Gyr (flat models) or 12 Gyr (open 
models), a baryon density parameter Qb ~ 0.0125/i^'^, and 
no influence of gravity waves or mild early reionization on 
the DMR anisotropies. With these assumptions, there is only 
a rather narrow range of DMR-normalized CDM models 
that produce the observed abundance of rich galaxy clus- 
ters. Flat models {Qo + Ao = 1) require Qo = 0.25 — 0.4, 
whereas open models require flo = 0.4 — 0.5. Most plausible 
departures from our assumptions tend to drive these esti- 
mates to higher values of Qo- For example, lowering Hq in 
order to increase to reduces the value of as implied by CMB 
anisotropies but has little effect on the estimate of as from 
the cluster abundance. Thus a lower Ho requires a higher 
fio, and fio = 1 is acceptable if /i « 0.25. Increasing the 
baryon fraction slightly above our standard value has an ef- 
fect in the same direction, but the impact is substantially 



weaker than the impact of lowering Hq unless Qb becomes 
a large fraction of Qq. Lowering the baryon fraction would 
go in the direction of allowing lower fio, but here the effect 
is quite small because our acceptable models already have 
a low ratio of Ob /Ho- Our adopted DMR normalizations, 
based on the COBE-DMR two-year data, tend to be higher 
(by roughly 10% or +l-a) than the central values estimated 
from the 4- year data (see discussion in Section 2.1). Low- 
ering the normalizations would reduce the value of as for a 
given flo, so again a higher flo would be required to match 
the observed cluster abundance. 

Generalizations of the simplest inflation models most 
easily lead to redder primordial power spectra, e.g., to an ef- 
fective index n < 1 in flat models. Lowering n reduces the as 
inferred from DMR-normalization both because it changes 
the extrapolation of fluctuation amplitudes to cluster scales 
and because, at least in power-law inflation models, such 
tilted spectra are accompanied by the production of gravity 
waves. Any change towards a redder inflation power spec- 
trum thus raises the lower limit on flo. Reducing n to about 
0.8 is a simple way to make flo = 1 compatible with the 
DMR data and the observed cluster abundance, though any 
ilo = 1 model remains difficult to reconcile with the high 
baryon fraction in rich clusters and with current observa- 
tional estimates of to and Hq. 

While our adopted range of median cluster masses is 
fairly conservative, we cannot exclude the possibility that 
this estimate may be systematically too low or too high. 
The gravitational lensing data suggest that, if anything, the 
former is more likely. In this case the observed cluster abun- 
dance would require a higher value of ug and hence a higher 
value of f2o for consistency with the DMR data. 

We conclude that, under the assumptions discussed 
above, Oo > 0.25 is required to produce enough mas- 
sive galaxy clusters in spatially-flat, DMR-normalized CDM 
models. In open models, the value of as implied by the DMR 
data drops so fast with decreasing fio that the lower limit 
rises to f2o ^ 0.4. The statistical significance of these limits 
is difficult to quantify because the errors are dominated by 
the systematic uncertainties just discussed. Our conclusions 
agree well with previous N-body work (Frenk et al. 1990; 
Bahcall & Gen 1993; White et al. 1993a; Eke et al. 1996b) 
and with semianalytic results based on the Press-Schechter 
model (White et al. 1993a; Eke et al. 1996b; Kitayama & 
Suto 1996; Viana & Liddle 1996). However, we are able to 
set a more stringent lower limit on f2o because the use of 
large N-body simulations removes any uncertainty regarding 
the relation between predicted and observed cluster masses. 
Related constraints on cosmological parameters follow from 
considering the evolution of the abundance of clusters (Eke 
et al. 1996b; Oukbir & Blanchard 1997). Current cluster evo- 
lution data tentatively indicate f2o ~ 0.4 — 0.5, as ~ 0.75 
(Carlberg et al. 1997a; Henry 1997 private communication), 
consistent with the constraints derived here for open models. 

There are possible changes to our assumptions or pa- 
rameter choices that would weaken these lower bounds on 
ilo. These include raising Hq (and consequently lowering to) 
or going to the high end of the error range for the DMR- 
normalization. Taking these and other uncertainties into ac- 
count, Gorski et al. (1996b) conclude that open- bubble in- 
flation models with fio as low as 0.3 could be consistent with 
current observations. Lower Q,q is also allowed if the inflation 
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model produces a bluer, "anti-tilted" primordial power spec- 
trum, though in some cases the effect of such a change on the 
DMR-normalized ag would be partly cancelled by the con- 
tribution of gravity waves to the DMR anisotropies. Open 
CDM models with n > 1 (on scales well below the Hubble 
length) have been advocated by Bunn & White (1997) and 
White & Silk (1996). In these studies the anti-tilted spec- 
trum is imposed ad hoc; it has not yet been shown that such 
a spectrum can be produced by quantum fluctuations in an 
open-universe inflation model. + 

Matching the CO-B S-DMR anisotropy amplitude and 
the abundance of rich clusters does not, by itself, guaran- 
tee a successful cosmological model. Of the many criteria 
available for discriminating amongst cosmological models, 
the observed pattern of galaxy clustering remains one of 
the most widely used. However, whereas the CMB and clus- 
ter abundance tests apply directly to the distribution of 
dark matter, the galaxy clustering test requires a model for 
the galaxy distribution. Unfortunately, our understanding of 
galaxy formation is still so primitive that only highly ide- 
alized mathematical models such as the "high peak" model 
of biased galaxy formation (Davis et al. 1985; Bardeen et 
al. 1986) have been explored in detail. Cosmological gasdy- 
namic simulations do not yet have sufficient resolution to fol- 
low galaxy formation except in relatively small volumes (Cen 
& Ostriker 1993; Katz, Hernquist, & Weinberg 1992; Sum- 
mers, Davis, & Evrard 1995; Frenk et al. 1996). Our dark 
matter simulations do not address the issue of galaxy for- 
mation, but they do at least provide an accurate description 
of the statistical properties of the dark matter distribution 
on scales ~ 0.2 — 35/i~^ Mpc. Comparing this mass distribu- 
tion with the observed galaxy distribution reveals the kind 
of biases that must be present in the galaxy distribution for 
a particular cosmological model to be acceptable. 

We find that all DMR-normalized CDM models that 
successfully reproduce the cluster abundance require that 
the galaxy distribution be biased relative to the mass in a 
non-trivial fashion. Except in the fio = 0.4, Ao = case, 
the acceptable low-fio models, with or without a cosmo- 
logical constant, produce mass distributions that are sig- 
nificantly more strongly clustered than the observed galaxy 
distribution on small and intermediate scales. These models 
thus require the galaxy distribution to be antibiased rela- 
tive to the mass on these scales. On scales ^ 10ft~^ Mpc, 
the mass correlations in these models become comparable 
to the galaxy correlations. In the f2o = 0.4, Ao = case, the 
predicted mass correlations are only slightly stronger than 
the observed galaxy correlations on small scales (and, within 
the uncertainties, could even be weaker), but they fall well 
below the observed ^(r) on scales larger than ~ 10fe~^ Mpc. 
This model thus requires (mild) antibias on small scales and 
positive bias on large scales. Finally, the tilted Oq = 1 model 
requires a strong positive bias everywhere, and the degree 
of bias depends on scale. Beyond ~ 10/i^^ Mpc, it must rise 



■f A recent study (Garcia-Bellido & Linde 1997) indicates that 
quantum fluctuations in an open-bubble inflation model with 
three scalar fields might produce an anti-tilted generalization of 
the usual single-scalar-field open-bubble inflation model primor- 
dial power spectrum (Ratra & Peebles 1994,1995; Bucher et al. 
1995; Yamamoto ct al. 1995). 



steeply with scale to compensate for the rapidly declining 
large-scale clustering strength. 

Little is known about the physical plausibility of a 
strong antibias in the galaxy distribution on small and inter- 
mediate scales, which would be required for almost all our 
low-f2o models to match observations. Galaxy mergers might 
produce this kind of effect, and preliminary results from 
a large programme of cosmological gasdynamic simulations 
suggest that galaxies may indeed be antibiased in low f2o 
models, although at a lower level than seems required by ob- 
servations (Jenkins et al. 1996), except possibly for f2o ~ 0.4 
open models. On empirical grounds, however, an antibiased 
galaxy distribution in models with f2o > 0.2 seems diffi- 
cult to reconcile with virial analyses of galaxy clusters. If 
the mass-to- light ratio of clusters equals the universal value, 
then these imply that f2o = 0.19, with a 2-a upper limit 
Q.0 < 0.33 (Carlberg, Yee & EUingson 1997b; Carlberg et al. 
1997c). However, if galaxies are over-represented in clusters 
by a factor B, then the data imply Q.q — 0.195. If galaxies 
were antibiased on small scales, then we would expect them 
to be under-represented in clusters, implying B < 1, and 
therefore Hq < 0.19 (or < 0.33 at 2-a). 

Given these considerations, our results lend support to 
earlier suggestions (Efstathiou et al. 1990; Stompor et al. 
1995; Ostriker & Steinhardt 1995; Ratra et al. 1997; Ganga 
et al. 1996c; Klypin et al. 1996; Maddox et al. 1996; Peacock 
1996) that low-density, flat-universe CDM models can more 
easily flt the observations if they have a mildly tilted primor- 
dial mass power spectrum. Flat-A models with n ~ 0.9 could 
still have reasonable values of fio and to for a given h, and 
they would require less small-scale antibiasing in the galaxy 
distribution. It is possible but less clear that open CDM 
models would benefit from an anti-tilted primordial spec- 
trum, for although this would allow a lower f2o and higher to 
for a given h, it would also lead to stronger small-scale mass 
clustering and hence a requirement of substantial antibias in 
the galaxy distribution on these scales. A bluer primordial 
power spectrum could also bring the small-scale (multipole 
(. > 300) CMB anisotropies in open models into conflict with 
upper limits (Tucker et al. 1993, Ratra et al. 1997, Church et 
al. 1996, Ganga et al. 1996a), and it might result in a steeper 
CMB anisotropy spectrum on intermediate scales {£ ~ 200) 
than is indicated by current data (Ganga et al. 1996c). Both 
of these effects, however, could be ameliorated by mild early 
reionization. 

To illustrate the region of parameter space allowed by 
current observational uncertainties, consider an flo = 0.35, 
open-bubble inflation model with to = 12 Gyr (as advocated 
by Alcock et al. 1996) and Qsh^ = 0.007 (as advocated by 
Songaila et al. 1994; Carswell et al. 1994; Rugers & Hogan 
1996a,b). At the 2-a upper limit, the DMR-normalized ag 
for this model is 0.69 (Gorski et al. 1996b), consistent, also 
within 2-a, with the Eke et al. (1996b) cluster abundance 
requirement and with the constraint from cluster mass-to- 
light ratios that B > 1.1 {2-a, Carlberg et al. 1997b,c). 
The slightly smaller fio compared to the 04 simulation il- 
lustrated in Fig. 5 would lead to a slightly lower F and, on 
large scales, this change would compensate for the reduc- 
tion in clustering strength caused by reducing erg from 0.75 
to 0.69. An open model with these parameters might thus be 
consistent with the DMR anisotropies, the evolution of the 
abundance of clusters, the APM clustering data, and cluster 
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mass-to-light ratios if it had a weak positive bias on smaU 
scales and a somewhat stronger positive bias on larger scales. 
Agreement at this level could probably also be achieved by 
a tihed, flat- A CDM model (Ostriker & Steinhardt 1995). 

A great variety of observations currently underway or 
planned for the near future will significantly sharpen the 
constraints that we have discussed in this paper. Analyses 
of existing small-angle CMB anisotropy data, based either 
on comparisons of observations from many different experi- 
ments with model predictions for a range of parameter val- 
ues (e.g., Ratra et al. 1997; Ganga et al. 1996c; Hancock 
et al. 1996; Lineweaver et al. 1997; Rocha et al. 1997), or 
on maximum likelihood analyses of CMB data sets that 
directly use model CMB anisotropy spectra (e.g., Ganga 
et al. 1996a,b; Bond & Jafi'e 1996), can be used to con- 
strain parameters, although they do not yet provide strong 
constraints in the most interesting regions of parameter 
space. Higher precision measurements of CMB anisotropics 
will become available in the next few years, from numer- 
ous ground- and balloon-based experiments and eventually 
from the MAP and COBRAS/SAMBA space observatories. 
These measurements can, in principle, constrain cosmologi- 
cal parameters to unprecedented accuracy. At the same time, 
our knowledge of galaxy clustering will improve significantly 
with the next generation of redshift surveys, the SDSS and 
2dF surveys (Gunn & Weinberg 1995; CoUess 1996). (The 
simulations described here will be used to generate pub- 
licly available mock catalogues with the geometry and se- 
lection criteria appropriate to these surveys.) The ability 
to model galaxy formation using cosmological gasdynamic 
simulations is improving rapidly, and it may well soon be 
possible to establish whether the apparently complex rela- 
tion between galaxies and mass implied by our analysis is a 
natural outcome of hierarchical clustering in a low-fio uni- 
verse. 
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